Axially symmetric Einstein-Straus models 

Marc Mars* 

School of Mathematical Sciences, Queen Mary and Westfield College, 
Mile End Road, London El 4NS, United Kingdom. 

The existence of static and axially symmetric regions in a Friedman-Lemaitre (FL) 
cosmology is investigated under the only assumption that the cosmic time and the 
static time match properly on the boundary hypersurface. It turns out that the most 
general form for the static region is a two-sphere with arbitrarily changing radius which 
moves along the axis of symmetry in a determined way. The geometry of the interior 
region is completely determined in terms of background objects. When any of the most 
widely used energy-momentum contents for the interior region is imposed, both the 
interior geometry and the shape of the static region must become exactly spherically 
symmetric. This shows that the Einstein-Straus model, which is the generally accepted 
answer for the null influence of the cosmic expansion on the local physics, is not a robust 
model and it is rather an exceptional and isolated situation. Hence, its suitability for 
solving the interplay between cosmic expansion and local physics is doubtful and more 
adequate models should be investigated. 
PACS numbers: 04.20.Jb, 04.20. Cv, 98.80.-k, 04.40.Nr 



*Also at Laboratori de Fi'sica Matematica, Societat Catalana de Ffsica, IEC, Barcelona. 



1 



I. INTRODUCTION 



The influence of the large scale geometry of the Universe on the local physics around astrophysical 
objects is a very old and fundamental question in gravitational physics. The fact that the geometry 
of the Universe as a whole can be approximated very accurately by an expanding Friedman-Lemaitre 
(FL) model leads naturally to consider whether the cosmological expansion has any effect on the local 
physics at astrophysical scales. This type of question was probably first addressed by McVittie M who 
used a spherically symmetric spacetime supposedly describing a point particle immersed in a Friedman- 
Lemaitre background (see, however, g]). This metric was later analyzed in more detail by Jarnefeld 
H . Both authors concluded that the effect of the cosmological expansion on the planetary orbits is 
very small and with no measurable consequences. More recently, Gautreau Q analyzed the problem 
by using an inhomogeneous and spherically symmetric model in the so-called curvature coordinates 
|^| (see H for a review of these results). However, the generally accepted answer to this problem was 
given by Einstein and Straus in Q| where a model describing a massive point particle surrounded by 
a spherically symmetric vacuum region embedded in a dust FL cosmology was presented (the now so- 
called Einstein-Straus model). The boundary of the vacuum region is a two-sphere comoving with the 
cosmological flow. Thus, Einstein and Straus showed the possibility of having static regions embedded 
in an expanding cosmological background thus implying that the cosmic expansion does not influence 
the local physics around the massive object. Despite its clear implications and simplicity, Einstein 
and Straus' construction suffers from several important limitations. First, the interior geometry of 
the static region is exactly spherically symmetric and its boundary is a two-sphere comoving with the 
cosmological flow. This a very strong assumption, specially regarding the shape of the static cavity. 
Second, given the dust background metric, the radius of the vacuum region surrounding the compact 
object is uniquely fixed by its mass. This implies that some situations of astrophysical interest are not 
satisfactorily described by the Einstein-Straus model. Furthermore, Krasihski has argued [|| that the 
model is unstable, meaning that a matching between the vacuum cavity and the FL background with 
a slightly bigger radius than the Einstein-Straus value makes the static region expand with respect to 
the cosmological flow (and vice versa for a smaller radius). 

Hence, more general models should be analyzed in order to reach a definitive conclusion on whether 
the influence of the cosmological expansion on the local physics is indeed vanishing and, in that case, 
on the suitability of the Einstein and Straus' construction to provide a completely satisfactory expla- 
nation. In order to test the importance of the spherically symmetric assumption, other geometries for 
the interior region (especially the shape of the boundary) should be considered. Since the local geom- 
etry is generally believed to be static (which is the translation of the null influence of the background 
expansion), the problem can be formulated as to whether non-spherical static cavities can be embed- 
ded in a FL cosmology. In a very recent paper, Senovilla and Vera ||] proved the surprising fact that 
embedding a cylindrically symmetric static region in an expanding FL cosmology is always impossible 
(irrespective of the matter contents inside the cavity). This result supports the view that spherical 
symmetry may indeed be an essential ingredient for the existence of the static region. In order to 
clarify this question, we address the existence of static cavities in a FL background in a rather general 
situation. The essential simplifying assumption we make (otherwise the problem becomes formidable) 
is that the interior static region is axially symmetric, both in shape and in spacetime geometry (no 
restrictions on the interior energy- momentum contents are imposed). This assumption is physically 
very reasonable because axial symmetry is likely to be present in many realistic situations. The second 
assumption we make is that the two canonical time functions in each region (the cosmic time in the 
FL background and the static time in the interior region) match properly across the boundary of the 
static region, so that the global spacetime has a well-defined and natural time function. This model is 
general enough to include both spherical symmetry and cylindrical symmetry as particular cases and, 
hence, its analysis will clarify to what extent the extreme situations in which the matching is possible 
(spherically symmetry) and never feasible (cylindrical symmetry) are generic. 

Considering static regions embedded in an expanding background requires the use of the theory of 
matching between spacetimes. This consists in imposing the matching conditions which can be viewed 
as a set of partial differential equations to be solved. When the symmetry of the matching problem 
is high enough (i.e. when the two spacetimes to be matched possess at least two isometries which are 
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preserved by the matching hypersurface) , the equations become ordinary differential equations and 
the problem simplifies considerably. Most of the explicit matchings which have been performed in the 
literature allow for this type of reduction. In particular, this is true for the Einstein-Straus model and 
the cylindrically symmetric case described above. In our case, however, the problem is technically 
more complex because the global spacetime possesses only one isometry, namely axial symmetry. In 
order to solve it, we use a set of constraint matching conditions, which are necessary consequences 
of the matching conditions on any spacelike two-dimensional surface embedded into the matching 
hypersurface. They appear naturally when developing a 2 + 1 geometrical decomposition of the 
matching conditions into a set of constraints on initial spacelike two-surfaces and a set of evolution 
equations determining how the two-surface must evolve in order to generate the three-dimensional 
matching hypersurface. The full decomposition has not been completed yet and will be reported 
elsewhere. For the purposes of this paper it will suffice exploiting the set of constraint matching 
conditions. 

The paper is organised as follows. In section 2, the matching theory between spacetimes is briefly 
reviewed and the constraint matching conditions on spacclikc two-dimensional surfaces are presented. 
In section 3, the matching problem between a static and axially symmetric spacetime and a Friedman- 
Lemaitre expanding cosmology is solved. Since the problem is technically non-trivial, the procedure 
is described in some detail. First, the constraint matching conditions are imposed. They restrict the 
form of the matching hypersurface severely. It turns out that the most general matching hypersurface 
consists of a collection of two-spheres, two-planes or hyperbolic planes centred on the axis of symmetry 
and moving arbitrarily along this axis. For closed FL models only two-spheres are possible; for flat 
FL models both two-spheres and two-planes are allowed while for the open case all three cases are 
possible. We restrict the analysis to the two-sphere case because this is the only one describing a 
spatially compact region embedded in a FL cosmology. Then, the remaining set of matching conditions 
are imposed. They determine the interior geometry of the static region in terms of the background 
FL metric and the function r(t) describing how the radius of the two-spheres generating the matching 
hypersurface varies with time. The motion of the set of two-spheres along the axis of symmetry 
is also determined by the background metric and r(t). In section 4, a detailed study of the static 
and axially symmetric interior geometry is performed and a number of important conclusions are 
obtained. In particular, the interior metric is shown to be nearly spherically symmetric with both 
its energy-momentum tensor and Petrov type of the same type as in the exact spherically symmetric 
case. Furthermore, analyzing the most widely used energy-momentum tensors for the interior region 
(vacuum, A-term, perfect fluid and electrovacuum) we show that the interior geometry must be exactly 
spherically symmetric in all those cases and the boundary of the cavity must be generated by a 
two-sphere with its centre being at rest with respect to the cosmological flow. Hence, we conclude 
that the Einstein-Straus model is essentially unique (when axial symmetry is imposed) and that 
the spherical symmetry assumption is indeed crucial for the existence of static cavities inside a FL 
cosmology. Consequently, obtaining a robust model determining whether or not the cosmological 
expansion has any influence on the local physics is still an open problem which deserves further 
analysis and investigation. 



II. JUNCTION OF SPACETIMES AND CONSTRAINT MATCHING CONDITIONS. 

The matching (see j^] and references therein) between two spacetimes with boundary, (V + , g + ) and 
(V~,g~) with corresponding boundaries fl + and £l~ , consists in constructing a spacetime manifold 
(V, g) by identifying the points, (and the tangent spaces) on the two boundaries, which are therefore 
called matching hyper surf aces. The point-to-point identification requires the existence of a diffeomor- 
phism between fi + and Q - . Therefore, they can be also considered as diffeomorphic to an abstract 
three-dimensional manifold Q so that there exist two embeddings 

<J>+ : n — ► V + , $~ : fi — > V~, 

satisfying $+ (0) = 0+ and $~ (fi) = Q~. Using local coordinates for Q (i,j, ■ ■ ■ — 1, 2, 3), {x+ } 
for V + (a,0, ■ ■ ■ = 0, 1, 2, 3) and {x" } for V~ the two embeddings take a local form 
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where (£) and x°_ (£) are functions of £*. As an important theorem by Clarke and Dray jlO| shows, 
the construction of a spacetime V = V~ U with continuous metric g can be performed if and only 
if the two first fundamental forms <& + (g + ) and <f> _ (g) on f2 coincide ($* denotes the pull-back of 
the map $). In local coordinates, these conditions read 
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and will be called preliminary matching conditions. When ([!]) is fulfilled, the Riemann tensor in (V, g) 
is well-defined in a distributional sense. In general, this distribution has a Dirac delta part with 
support on the matching hypersurface, which can be used, when appropriate, to describe thin shells 
of matter and/or pure gravitational field. Whenever these physical objects are absent, the Dirac delta 
contribution must vanish and the second set of matching conditions must be imposed. First, two 
transverse vector fields, 1+ on Q + and I- on Q~ , (called rigging vector fields) must be chosen. As 
they prescribe how the tangent spaces are to be identified, their choice must be compatible with the 
continuity of the metric in V . This demands 
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The two riggings and i_ must have different relative orientations with respect to f2 + and f2 , 
respectively (i.e. either 1+ points outside V + everywhere and I— points inside V~ everywhere or 
vice versa). Once the riggings are chosen, the matching conditions demand that the pull-back of the 
covariant derivative of the rigging one-form (obtained by lowering the index to I ) coincides when 
calculated using any of both embeddings. In local coordinates, 
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where V + (V ) denotes the covariant derivative and (^uv) are the corresponding Christoffcl 
symbols in V + {V~). 

In many physically interesting problems the degree of symmetry is high enough so that the con- 
ditions (0) and (^) become a system of ordinary differential equations, which simplifies the problem 
considerably. In general, however, this reduction does not happen and the partial differential system 
must be considered. As we shall see below, some of these cases can be handled more easily by using 
a 2 + 1 decomposition of the matching conditions. This consists in obtaining a set of geometrical 
conditions on an initial two-dimensional surface and evolution equations determining how this surface 
must evolve in order to generate the matching hypersurface. This theory will be reported elsewhere 
when fully developed. For the aim of this paper, only the constraint matching conditions will be 
necessary. In order to describe them, some notation must be introduced. 

The main assumption for the 2 + 1 decomposition is that the hypersurface fi can be foliated by a set 
of two-dimensional surfaces, which will be denoted by {S T } (r is a parameter identifying the element 
of the foliation and plays the role of an evolution parameter). Each two-surface in the foliation can 
be immersed into f2 by using the natural inclusion 
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The surface S T can then be embedded into V + and V ~ , respectively, by the maps 

$+ = $+oj r , $~=$~oi T (4) 

which define the two image surfaces £ + = $ + (SV) and S~ = (SV). We will further assume that 
both S + and S~ are spacelike everywhere. Let us consider a point i£S T and define the vector space 

(T x St) x ^{ne Tl t[x) V + ■ n (T X S T )) = o} , 

where d$ + is the differential application of the embedding $+, and T y M, T*M denote, respectively, 
the tangent and cotangent spaces of a manifold M at any point y e M. The vector space (T x S+)j_ is 
simply the set of normal one- forms to the two-surface S + in y + . Since S + is assumed to be spacelike, 
(T x S^)j_ is two-dimensional and timelike. The second fundamental form of S + along the direction n 
is defined, as usual, as 

(n) = $+* (V+n) . 

where n is a one-form field orthogonal to S + everywhere (the set of all these objects will be denoted 
by (TS + )j_). Obviously, similar definitions hold by interchanging plus and minus everywhere. 

The constraint matching conditions on S T are essentially the restrictions of the full set of matching 
conditions into the two-surface S T . They can also be splitted into two different sets, one involving first 
fundamental forms and another involving second fundamental forms. The first one is an immediate 
consequence of (Q) and the definition (Q) and demands, quite naturally, the equality of the two first 
fundamental forms on S T . Explicitly, 

$r (.9 + ) = *v* (.9-) • (5) 

The second set of constraint matching conditions is somewhat more involved, although its geometrical 
contents can be clearly understood. The matching of the two spacetimes, requires, as mentioned above, 
the identification of the tangent spaces at the points on the matching hypersurfaces. Furthermore, 
this identification must be done so that the resulting metric tensor is continuous. When working with 
the foliation {5V}, part of the job is accomplished by identifying the tangent planes to the two surfaces 
S + and S~. The remaining identification requires the existence of a linear and isometric application 

(T X S+) X — (T x S-) ± , (6) 

which completes the identification of the tangent spaces at the points $ + (x) and <£>~(x). ft is not 
difficult to show that the preliminary and second set of matching conditions can be combined to imply 
the following condition on S r 

K+ r (n) = X- (£ r (n)) , Vn e (TS+) ± , (7) 

where the field £ T (n) is defined, naturally, as £ r (ro)|$-( x ) = (n\^+^ x -^, Vx e S T . The existence of 
an isometric application (^) satisfying (|7]) constitutes the second set of constraint matching conditions. 
Due to the linearity of and the ^"-linearity of Kg and (i.e. (/n) = fK^ (n) for any 

scalar function /), it suffices to impose (Q) for two one-form fields ni,ri2 G (TS+) X being linearly 
independent at each point on S + . Let us now write down the local form of the constraint matching 
conditions (||) and (0). Introducing local coordinates {( A } (A,B,--- — 1,2) in the two-surface S T , 
the two embeddings read, locally, 

C A ^< = * +a (C), K - ^ A ^x1=x- a (0, 

where x+ a (£) and x~ a (£) are functions of C A - The condition (||) becomes, obviously, 
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Regarding <fj\), we take arbitrary bases, n4 of (T x S+) ± and ni 
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the scalar quantities 
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of (T X S T and define 



The linear isometric application £^ takes the form 
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where ( T ^(i) are a set of scalar functions defined on S T which must satisfy 

^r A B {x)^ C D {x) 1 BD {x)^ 1 i C { X ) 

to ensure that is an isometry. The condition (Q) can now be rewritten as 
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In the next section, the conditions (ra) and (^) will be exploited to restrict the form of the matching 
hypersurface in the axially symmetric Einstein-Straus model. They will allow for a complete resolution 
of the matching problem we are considering. 



III. AXIALLY SYMMETRIC EINSTEIN-STRAUS MODEL. 



Let us now study the main object of this paper, namely the existence of static, axially symmetric 
regions embedded in a Friedman-Lemaitre universe. We choose standard spherical coordinates in the 
FL cosmology so that the line-element takes the usual form 

ds 2 FL = -dt 2 + a 2 (t) [dr 2 + £ 2 (r, e) (d9 2 + sin 2 9d(j) 2 )] , (10) 

where a(t) is the scale factor of the cosmological model (we do not restrict its form whatsoever a 
priori) and £ (r, e) is the standard function 

{sinhr if e = — 1 
r if e = . (11) 
sinr if e = 1 

Regarding the static, axially symmetric spacetime, we will make the usual assumption that it is 
orthogonally transitive (i.e. the two-planes orthogonal to the isometry group orbits are surface- 
forming). As an important theorem by Kundt and Trumper |llf shows (generalizing Papapetrou's 
theorem for vacuum fl2|| ), this condition is automatically satisfied for a wide class of energy- momentum 
tensors including vacuum, A-term (i.e. cosmological constant), perfect fluids and electrovacuum. Thus, 
the orthogonally transitive assumption does not represent a severe restriction on the kind of regions 
we are considering. Under this assumption, there exist coordinates in which the line-element in the 
static region reads 

ds 2 st = -F 2 (x A ) dT 2 + g AB (x c ) dx A dx B + W 2 (x A ) dj) 2 , (12) 

so that &t is the static and is the axial Killing vector. The coordinate transformations keeping the 
structure of (|l~2| ) unchanged are 

T — > a\T + ctQ, x A -> x A (x A ) , 0->0 + c, (13) 
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where ot\ ^ 0, «o & n d c are arbitrary constants and x A (x A ) are arbitrary functions with non- vanishing 
Jacobian. 

Both FL and static spacetimes admit canonical foliations by spacelike hypersurfaces. In the FL 
cosmology, this foliation is defined by the homogeneity and isotropy hypersurfaces t = const. In 
the static region, the foliation is given by the hypersurfaces orthogonal to the static Killing vector, 
T = const. Since the physical problem we treat is the existence of static regions living in a cosmo- 
logical background, we can assume that the matching hypersurfaces are not tangent to any of the 
hypersurfaces t = const, in the FL spacetime or T = const, in the static spacetime. This case would 
rather correspond to a phase transition between a static spacetime and an expanding FL cosmol- 
ogy, which is a completely different situation. Consequently, the foliation defined by the cosmic time 
induces a canonical foliation in the matching hypersurface fl FL , defined by S FL — Sl FL n {t = r}. 
Similarly, the matching hypersurface in the static spacetime f2 s * (the labels FL and St will replace 
the scripts plus and minus used in the previous section) can be canonically foliated by the two-surfaces 
S^f — fl St n {T — T}. Hence, we have that each region in the glued spacetime has a natural and 
canonical time function, the cosmic time in the FL background and the static time in the interior 
region. It is natural from the physical point of view to demand that these two canonical times coin- 
cide on the matching hypersurface. Here, we make the second fundamental assumption in this paper, 
namely, we demand that the global Einstein-Straus spacetime possesses also a canonical time function 
which coincides with the cosmic time in the FL background and with the static time in the interior 
region^]. If this assumption were violated, then spacetime events on the boundary of the static region 
being equally old with respect to the cosmic time, would correspond to different instants of time with 
respect to the static time, which is an undesirable property for the model we are constructing. Notice 
that this condition we are imposing is satisfied in both the spherically symmetric Einstein-Straus 
model and in the static cylindrically symmetric model in raj. Therefore, the assumption is not only 
physically plausible but it is also general enough so that a bridge between the two extreme cases is 
constructed. Hence, analyzing this case will certainly determine whether the spherically symmetric 
Einstein-Straus model is exceptional or not. It must be emphasized, however, that we are indeed mak- 
ing a simplifying assumption which is not mathematically necessary a priori. It would be interesting 
analyzing the problem in full generality, although this is probably a much more difficult question. 

The condition we are imposing can be translated in mathematical terms by demanding that the 
two canonical foliations in the matching hypersurface coincide in the glued manifold. In other words, 
we assume there exists a foliation {S T } in the abstract matching hypersurface f2 such that the two 
embeddings defining the matching hypersurfaces 

$ FL : Q — > V FL , $ st : tt — ► V St , 

satisfy <I> FL (S T ) = S FL and <fr s * (S T ) = S^ T \ (the geometrical coincidence of the foliation allows for 
a reparametrisation T(r)). 

Regarding the condition that the static region is axially symmetric, this implies not only that the 
interior metric possesses an axial Killing vector, but also that the two matching hypersurfaces preserve 
this symmetry. More precisely, they must be invariant under the action of the axial symmetry in each 
spacetime. Since the two spacelike foliations in V FL and V St are also axially symmetric, it follows 
that each two-surface S T admits an axial Killing vector field. Choosing coordinates {A, (p} on S T 
adapted to this axial symmetry (i.e. such that the axial Killing reads d v ), the most general form for 
the two embeddings of S T into V FL and V St can be easily seen to be 



x If the interior region possesses a high enough degree of symmetry (e.g. Minkowski spacetime) the static 
time may be not unique. However, since we did not specify the form of the interior metric ( |l2| ) a priori, we 
are not choosing between any of them (whenever there is more than one) and therefore they are all treated 
on the same footing. In this particular case, the assumption we are making is that there exists at least one 
static time which matches appropriately with the cosmic time on the matching hypersurface. As we shall see 
below, the matching problem will fix the interior geometry completely and, therefore, it will choose which of 
the static times (if there is more than one) is the appropriate one to match with the exterior FL background. 
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$f : {A, ^} — {T = T(r), ^ = a£ (A) ,0 = ^} 



(15) 



where r r (A), # T (A), x A (A)and T(r) are unknown functions. Imposing the first constraint matching 
conditions (^|) gives the two relations 
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where the dot means derivative with respect to A. Regarding the second set of conditions (0), we must 
first identify the normal planes to the two-surface S T in each spacetime and then construct the most 
general linear isometry between them. The vector space (T^Sf*) . is spanned by the two mutually 
orthogonal one-forms 



dT 
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where cab is the totally antisymmetric symbol with ei2 = 1. Similarly, the vector space (T X S^ L ) ± is 
spanned by the two mutually orthogonal one-forms 



dt 



■S> FL (x) 



9 T dr — r T d6 



Thus, the most general linear and isometric application £^ : (T X S^ L ) 
the linear extension of 
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where (3 (A, 0) is an arbitrary function, 771 and r\2 are arbitrary signs and we have defined y 2 
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and G = det(<?As)- A consequence of the staticity of the interior region is 
= 0. Consequently, the second set of matching conditions @) becomes 

K s? =^^sinh/3K|* (/cU« (a)> 
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which implies 
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This equation involves only objects from the FL background metric. A trivial calculation shows that 
Kfj L ( K ^\^ FL (x)j an d ( K a \^ FL (x)j are both diagonal and therefore ( fl6| ) implies two conditions. 

Combining them so that the hyperbolic angle (3 disappears, we obtain the equation 
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where we have dropped a global factor a tt (comma means, as usual, partial derivative) which is non- 
zero because the cosmological model is assumed to be expanding (our aim is precisely analyzing the 
possible effect of the cosmological expansion in the local physics). (17) is an ordinary differential 
equation which defines the possible shapes of S^ L . It can be explicitly solved to give 
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where a\ and <jq are arbitrary integration constants (the particular solution 9 — tt/2 can be obtained 
from this expression by putting ctq = 1 and o\ — and following an appropriate limiting procedure). 
In order to identify this surface, let us analyze its intrinsic geometry. Using the defining relation ( |l8[ ) 
it is not difficult to obtain the first fundamental form in S^ L , which reads 
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Evaluating the curvature scalar Rsfl for this metric we find 
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which is constant throughout S^ L . Thus, this two-surface is maximally symmetric and hence p3[ 
is either a two-sphere, a two- plane or a hyperbolic plane depending on whether the sign of -R5F1, is 
positive, zero or negative, respectively. Due to fll9|), it follows that in a closed FL spacetime (e = 1) 
the matching two-surface S^ L must be a two-sphere, in a flat FL cosmology (e = 0) both two-spheres 
and two-planes are possible and in an open FL cosmology (e = —1) the three possibilities above are, 
in principle, allowed. The physical problem we are investigating is the existence of axially symmetric 
regions immersed in a Friedman-Lemaitre background. When the surfaces S^ L are two-planes or 
hyperbolic planes, the static regions are spatially unbounded (they extend to the spatial infinity in 
the FL cosmology) and the only possible construction for having a static region in a FL background 
would be performing a double matching, so that the static region is "sandwiched" between two FL 
spacetimes (i.e. the static region would have finite width but would be unbounded in the other two 
spatial directions). Even in that case, the static region cannot be prevented to disconnect the FL 
background. Thus, these two cases are not physically interesting for the problem we are considering 
in t his p aper. Consequently, we will restrict the analysis to the case in which {S^ L } are two-spheres 
(i.e (|19j ) is always positive). 

Let us now locate these two-spheres within the FL spacetime. To that end, we take advantage of the 
high degree of symmetry in the cosmological background. The uniparametric group of transformations 
generated by the Killing vector 



k = cos 9d r — sin 9da 



(20) 



is the most general translational isometry leaving the symmetry axis of invariant. Let us use the 
natural parametrization (which will be denoted by 7) induced by ( pp| ) in this group of transformations. 
It is not difficult to show that the transformation 7 = 70, where 70 is given by 
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transforms the surface (18) into the two-sphere r = ro, with ro given by 
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Thus, the matching hypersurface in the FL cosmology can be, at most, a two-sphere moving arbitrarily 
along the axis of symmetry and with a radius which can change freely as time evolves. At this 
stage, the result in || stating that no cylindrically symmetric Einstein-Straus models can exist, has 
been recovered and generalized. In those models the matching hypersurface consists of a collection 
of cylinders with arbitrarily varying radius. The results above show that any model in which the 
matching hypersurface at any instant of cosmic time is not a two-sphere, a two-plane or a hyperbolic 
plane is impossible. 

In order to complete the matching procedure, it is convenient to introduce new coordinate systems 
in the FL and static spacetimes so that the form of the matching hypersurfaces becomes simpler. To 
that end, we choose spherical coordinates on each hypersurface t — r in the FL spacetime so that 
the two-sphere S FL is centred at the origin. Completing these coordinates with the global time t, we 
obtain a well-defined system of coordinates. In order to obtain the line-element in the new coordinate 
system, we first change the original coordinates ( |l0| ) into cylindrical coordinates {t,p, z,(fi} by the 
standard transformation 



S (p, e) = £ (r, e) sin 



£ (z, e) £ (r, e) 



£, z (z,e) £ r (r,e) 



cosf 



(21) 



(t and 4> remain unchanged). In cylindrical coordinates, the Killing vector ( pOj ) reads simply k = d z 
and the corresponding uniparametric group of transformations takes the trivial form 



t -» £, 



z 2 + 7o(i), 



(22) 



where 70 is allowed to change in time because the centre of the two-sphere S^ L can move along the 
symmetry axis. Applying ( p2| ) and then transforming back into spherical coordinates by using (^l|), 
we obtain the linc-elcmcnt 



ds FL — 



-dt 2 + a 2 (t) 



(dr + f(t) cos 6dt) 2 + ( £ (r, e) dQ - /(i)S, r (r, e) sin 6dt) 

+S 2 (r, e) sin 2 



(23) 



where f(t) = 7o(£),< denotes the motion of the centre of the two-sphere along the symmetry axis and, 
for notational simplicity, the same symbols {£, r, 9, (f>} are used to name the new coordinates (which 
will be used from now on). The embedding defining the matching surface takes now the very simple 
form 



FL 



V 

{t = r, r = r(r), 6 = &, <p = <p} . 



(24) 



and the embedding <i> FL : S T — > V FL is obtained simply by fixing r = const, above. 

A convenient coordinate system in the static region can be constructed as follows. Due to the 
staticity of the spacetime, the integral lines of the static Killing 8t define a canonical diffeomorphism 
between different constant time hypersurfaces [Q. Thus, there exists a natural projection 



si 



n : V 

where 2 is the quotient space of V st defined by the orbits of the static Killing vector (it is, therefore, 
diffeomorphic to any of the hypersurfaces {T — const.}). Finding a coordinate system in an open 
subset U C S amounts to finding a coordinate system in the four-dimensional open set IT -1 (U) 
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because the static time T can be used as the fourth coordinate. Let us now consider the projection 
along II of the two-spheres {Sy*}. Since * is completely contained in the hypersurface {T = T(r)} it 
follows that fif* is isometric to II (5*^*) which, consequently, must be a two-sphere in S. Furthermore, 
S is axially symmetric (because of the axial symmetry of V St ) and each II (Sy ) is centred on the 
symmetry axis of S (this follows from the fact that {5*^*} are invariant under the action of the axial 
symmetry in V St ). Defining the subset of 5 



w = |J n (sf) , 

re J 

where J is the interval of variation of the parameter r, it is clear that either W is also a two-sphere 
or W has a non-empty interior. The first situation is marginal and it only happens when all the 
two-spheres {S^} are projected to each other by the integral curves of dr- Consequently, they must 
all have the same radius, which translates into E(r(r))a(r) = const, in the FL background. This 
fixes r(r) once the scale factor is given. Furthermore, a trivial calculation shows that the matching 
conditions ([!]) imply / = 0. Thus, the matching hypersurface in the FL spacetime is completely fixed 
and must be a two-sphere of constant geometrical radius with its centre being at rest with respect 
to the cosmological flow. Hence, the shape of the static cavity must be spherically symmetric in this 
case. 

In the generic situation, the set of matching two-spheres {5^* } are not all projected to each other by 

o o 

the integral curves of dr and, hence, W has a non-empty interior yy. At any point p <EW there exists 

o 

an open neighbourhood U p CW and a sufficiently small open subinterval I p C J such that the set of 
two-spheres {II (S^) , r £ Ip] foliates U p . In other words, for any point q £ U p there exists a single 

value of r q 6 I p such that q belongs to II [S^^j ■ A coordinate system in U p can now be constructed 

as follows. First, choose the natural angular coordinates {8, <f\ on each two-sphere LT (Sy *) such 
that 9 = and 9 = n correspond to the axial symmetry axis of S (this choice is possible because 
each two-sphere Sf* is centred with respect to the axis of symmetry in V st ). The third coordinate 
necessary to identify any point q £ U p C S can be chosen as the unique value T q £ I p identifying the 

o 

two-sphere passing through q. In this way, a coordinate system around any point p 6W is constructed. 

This coordinate system cannot, in general, cover the whole set W (because two different two-spheres 
with arbitrarily close values of r may intersect). Due to the purely local character of the matching 
conditions, this is not an obstacle to proceed with the matching. In order to simplify the matching 
equations and the final form of the interior static metric obtained below, it turns out to be convenient 
to use the symbol t to denote the coordinate identifying the two-sphere passing through any point 
q £ U p . This coordinate is spacelike by construction and it has nothing to do with the cosmic time in 
the FL background so that they should not be confused. The meaning of t should become clear from 
the context. The reason why the same symbol is used for such different objects is that both take the 
same value at every point on the matching hypersurface, which allows for an important simplification 
in the writing of the expressions below. 

The line-element for the static spacetime in the coordinate system just constructed takes the form 



ds st 



-F 2 (t, 9) dT 2 + Q 2 (t, 9) dt 2 + f 2 (t) 



d9 + H (t, 9) dt) + sin 2 9d$ 



(25) 



which is the most general metric for a static, axially symmetric and orthogonally transitive spacetime 
such that the two-surfaces T = const, t = const, are two-spheres. 

The embedding of the matching hypersurface into V St takes now the simple form 



st 



V 



si 



St 



{r,tf,yo} _> {t = T(t), t = r, = 0, = yj 



(26) 



and the corresponding embeddings : S T — * V are obtained by fixing r = const. 
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We are now in a favourable position to complete the matching problem. We already made use of the 
constraint conditions (||) and (|t]) in order to restrict the form of f2. However, not all their information 
has been extracted yet. The equality of the two first fundamental forms on S T gives now the single 
relation 



r=r(t) 



f(t) = a(f)£(r,e) 
The vector spaces (T X S? L ) ± and (T X S^) ± arc 

(T X S^ L ) ± = i^dt^FL^ , dr\^FL^ , (T x S^ t ) ± = (^dT^st^ , dt\ 9 st^ 

and the most general linear isometry between them is the linear extension of 

Cr( rf *l*^(x)) = Vi coshf3FdT + r/ 2 smh/3Qdt 



(27) 



/ sinh/3 



= //i F I - - / cos 6 cosh/3 ) 



(28) 



cosh/3 



/ cos 9 sinh /3 I eft 



*St(x) 



A straightforward calculation shows that (0) are equivalent to the three equations 



sinOHfi — cos 0H\ S = 0, 



772^ sinh j3 (Hf cos 9 — f ,< sin 6) — QT, 2 aa t t sin ( 
Ea,* cosh /3 + £ ,, sinh /3 



which can be trivially solved to give 

H(t,9) = L(t)smt 



Q(t,e) = 



Lf cos 9 



sinh/3 



-Ea, t 



El - S 2 a 2 



cosh/3 



r^/E2.-EV> 

Xj 7- 



r(t) 



El - S 2 C 



r=r(i) 



r(t) 



(29) 



(30) 



where L{t) is an arbitrary function. Thus, the 9 dependence of the two metric coefficients Q and H 
in (|25| ) is completely fixed. These expressions show that a necessary condition for the existence of the 
matching is 



eV, 



r(t) 



>0. 



(31) 



This inequality can be geometrically interpreted by noticing that the two-sphere S£ is non-trapped 
(see e.g. ]l5| for a definition) if and only if (|3l]) is fulfilled. Thus, it follows that a trapped (or 
marginally trapped) surface cannot be the boundary of an axially symmetric static region, which is 
an intuitively clear result. Using the expression ( |43| ) below for the cosmic energy-density p FL , the 
inequality (pin can be rewritten as 



E(r(f),e) < -I- 



a(t)Vp Fi («) 



(32) 



Therefore, the radius of the static region must be smaller than some maximum radius. In particular, 
this forbids situations in which the radius of the two-sphere S FL grows unboundedly until reaching 
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infinity (so that, in fact, it becomes a two-plane). Hence, if the matching hypersurface is a two-sphcrc 
at any instant of time, it must remain a two-sphere for ever (of course, this two-sphere can become a 
point and then disappear). 

Having exhausted the constraint matching conditions, we must now consider the more complicated 
set of conditions ([j]) and (||). The task is, however, much simpler now because of all the information 
we already know. The equality of the two first fundamental forms on provides 



£ f 

~e" 



/(*), 



FT, 



aa^S (r t + / cos 6) 



r=r(t) 



£ 2 r - a 2 £ 2 



(33) 



r—r(t) 



which fixes the 9 dependence of the remaining metric coefficient F in the static spacetime. Using (| 
( |29| ) and (33), the static line-element can be rewritten in the form 



ds st = 



[£ r + aa <t E(r it + fcos6)Y 



-a 2j 



cie 



■dT 2 - 



[a t £ + a£ r (r t + /cos6>)] ^ 2 



/ sin 9dt 



2v^2 • 2 

a Ij sin 



(34) 



r—r{t) 



which is completely given in terms of cosmological quantities, except for the function T(t) (which fixes 
the form of the matching hypersurface as seen from the static region) . In order to impose the second 
set of matching conditions we choose the rigging 



FL 



FL 



(f cos 6dt + dr 



a 

1st 



= r)iF sinh /3dT + r] 2 Q cosh fidt 



where the one- form Igt has been obtained from Ipi, by applying the linear isometry (p8[). Evaluating 
(Q) involves now a rather long calculation which was performed with the aid of the computer algebra 
program REDUCE Jl6| . Only the (r, r) component is not identically satisfied and it provides two 
conditions. Whenever / is not identically vanishing these two equations can be rewritten as 



Tuf = T t 



f,t + V 



a t t £ 2 r + a 4 ££, r r t (e + a 2 t ) 



(aa,t£r,i + E r ) a, t f.t + aa 2 t Y. tr f 3 + 



a t (E 2 r - a 2 E 2 ) 
[-aa 2 t (Er jtt + E, r r 2 t ) + 2 (a «S, r + ea^Er t )] / 



r=r(t) 



r—r(t) 



0. 



(35) 



(36) 



The case / = can be obtained from these equations by substituting / jt from ( pq ) into (|3^), then 
dropping a global factor / and finally putting / = 0. Only one equation survives in that particular 
case. 

The expressions ( j35[) and (|36| ) are differential equations for T(t) and fit). Thus, they determine 
the form of the matching hypersurface in the static spacetime and the motion of the centre of the 
two spheres {Sf L } along the axis of symmetry in the FL spacetime. Both equations can be explicitly 
integrated once the scale factor a(t) and the evolution r(t) of the radius of the two-spheres {S[ L } are 
given. In order to solve (|36|) two different cases must be considered. In the particular case 



aa (Er ( 



r(t) 



0. 



(37) 



the equation becomes algebraic and its solution is (the other solution / = is impossible due to (p3|)) 



aa*t 



a 2 a 2 t 



(P + P) 
2a 2 



FL 



13 



where the last equality follows from the standard expressions for the energy-density and pressure in 
the FL spacetime in terms of the scale factor (se e (|43|) below). 

In the generic case, (37) does not hold and (|36|) is & Bernouilli differential equation which can 
integrated U% in terms of quadratures as 



f fldt 



fit) = -j= 

Vl-2a 2 / f 2 e 2 J fldt dt 

where a is an arbitrary constant and f± (t) , ji (t) are defined as 



(38) 



hit) = 



Sr., 



2 (a tt E )T . + ea, t Er t ) 



a t (aa t Er t + £, r ) 



r=r(t) 



,h(t) = 



— aa^E^, 



aa.t 



Er t + £ r 



-=r(t) 



Regarding (|35|), this is a linear differential equation for T.t which is therefore trivial to integrate. The 
structure of its general solution is T(t) = aiT (t) + a$, where T (t) is any particular solution and ao, 
ot\ are integration constants. These constants are irrelevant as the form of the matching hypersurface 
remains unaffected if they change. They simply correspond to the freedom ( [l3| ) in defining the global 
static time in any static spacetime. Thus, ao and ct\ can be fixed arbitrarily without loss of generality. 

Thus, we conclude that the matching of a Friedman-Lemaitre spacetime with an axially symmetric 
static region allows for an arbitrary scale factor in the FL spacetime and for an arbitrary function 
r(t) describing the evolution of the radius of the matching two-spheres (as long as the inequality ( |32"| ) 
is fulfilled). Then, the matching conditions fix the motion of the centre of the two-spheres along the 
axis of symmetry through (|3^). They also fix uniquely the form of the matching hypersurface in the 
interior static region and the full geometry of the static spacetime (in an open region around the 
matching hypersurface), which is given by the line-element (34). In the next section, the geometry of 
the interior static geometry will be analyzed. 



IV. GEOMETRY OF THE STATIC SPACETIME. 

Let us first of all discuss the domain of applicability of the coordinate system {T, t, 6, cj)} in the 
static spacetime. As remarked in the previous section, these coordinates are well-defined at least in 
a sufficiently small open neighbourhood around any point in the interior of II~ 1 (n(£!' s *)). A priori, 

o 

there is no reason to expect this coordinate system to cover the whole open set II (w). However, this 
turns out to be the case and a single coordinate patch can be used to describe all the relevant interior 
region. Before showing this explicitly, it is worth pointing out that the form of the interior static 

metric at points not belonging to II (w) can not be determined at all by the matching conditions. 
This is because the matching procedure is purely local and there is no way of relating points outside 
this open set with the matching hypersurface. Determining the interior metric at those points requires 
further information like analyticity conditions or the Einstein field equations for some type of energy- 
momentum contents. Thus, all the discussion below concerning the geometry of the static region is 

o 

restricted to II -1 (yy)- 

In order to see explicitly that {T, t, 9, </>} is a valid coordinate system everywhere, we must show 
that no coordinate singularities are present in the line-element ( |34j ) or, equivalently, we must proof 
that any possible coordinate singularity is, in fact, an essential singularity. The metric ( |34| ) becomes 
singular at any point satisfying 

Y. <r + aa it Z (r it + f cos 6) =0 or a t E + a£ >r (r t + / cos 6) = 0. (39) 

The remaining possibilities, like T t = or / — > +oo, have clear implications for the matching 
hypersurface as seen from the FL spacetime (where the coordinates have a transparent meaning) and 
can be discarded as unphysical for the problem we are considering. In order to show that points 
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satisfying (|39| ) are essential singularities, let us write down the energy-momentum tensor of the static 
spacetimc. Using the orthonormal tetrad 



f aa. t £ (r tt + f cos ( 



-aT 



a jt £ + a£ r (r,i + / cos #) 



a£ d(9 



=r(t) 
£. 



dt 



£ 2 r - a 2 £ 2 



(40) 



-/sin 



a£ sin 



r—r(t) 



=r(t) 



r— r(t) 



to evaluate all tensor components, it turns out that the energy-momentum tensor is diagonal and 
with the structure (T a p) = diag (p,p r ,pT,pr), where p stands for the energy-density, p r is the radial 
pressure and pt the tangential pressure. We borrow the terms radial and tangential pressure from 
the spherically symmetric case even though the static metric does not possess this exact symmetry 
in general. However, the static metric is nearly spherically symmetric (its departure being governed 
by the single function /) since at every point we have a two-sphere passing through it. The structure 
of the energy-momentum tensor is exactly the same as in the exact spherically symmetric case. This 
resemblance goes even further because the evaluation of the Weyl spinor in the null tetrad associated 
canonically with (^) shows that 4"2 is the only non- vanishing component, exactly as in the spherically 
symmetric case. Thus, the static spacetime is algebraically special of Petrov type D (or conformally 
flat when * 2 =0). 

The explicit expressions for p and p r can be read from the expressions 



3a£ r (a 2 + e) (/ cos 8 + r, t ) + a, t £ (2aa,« + a 2 + e) 



a 2 [a t £ + a£ r (r >t + f cos 6)] 



a(p + p) FL (£ 2 r -a 2 £ 2 ) (/cos0 + r, t ) 



[a, t £ + a£ r (r, t + / cos 6)] [£ r + aa t £ (r t + / cos 9)] 



=r(t) 



r—r(t) 



(41) 
(42) 



where p L and p FL stand for the cosmic energy-density and pressure in terms of the scale factor 



P 



FL 



3 (a 2 +e) 



P 



FL 



-2aa 



(43) 



Regarding pt, its expression takes different forms depending on whether aa it Er jt 
the first case we obtain 



S, r ^ or not. In 



(p + pf t L a££ r + 2 (p + p) r " a£ 2 r (r, t + / cos 



FL 



p T 



{p + p) FL ^y. 



2 [a (t S + a£^ r (r >t + / cos 6)] [S jr + aa, t S (r t + f cos 6»)] 
aEE, r r,t (3a 2 t - 2aa,«) + a, t [E 2 r (3 + a 2 / 2 ) - a 2 (l + eE 2 ) r 2 - a 2 SE, r r, i( ] 



2 [a t £ + a£, r (r t + f cos 6»)] [£ r + aa, t £ (r, t + / cos 6>)] [£^ r + aa t £r t ] 



r— r(f ) 



and in the second case (aa ; t£r it + £ 
P 



0), the corresponding expression reads 



Pt 



4aa, t £/cos6' (a 2 £ 2 - £ 2 , + aa, t ££, r /cos 
+ 2 (p + p) FL £ r [2aa, t ££, r / cos 6 + £ 2 



^ aa :t £ 2 £ :I . (p + p) 
+ 4e - ao,«) - 2] 



t 



(44) 



Since the scalar quantities p, p r and pt have a clear physical meaning, any divergence in them 
corresponds to a true singularity in the static metric. It is now clear that at any point satisfying ( |39| ) 
we have divergences in p and/or p r and, therefore, we have a matter singularity rather than just a bad 
coordinate description. This type of singularities leads to physically unrealistic models and should, 
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therefore, be discarded (notice that these singularities are "touched" by the matching hypersurface and 
cannot be interpreted as singular "sources" for the static metric, like in the Schwarzschild spacctime) 
The expressions for p, p r and pt above allow also for obtaining a number of interesting conclusions. 
When the FL spacetime is de Sitter (or anti-de Sitter) then (p + p) FL = and the energy- momentum 
tensor in the static spacetime becomes T a p = Ag^ (with the same value for A as the cosmic back- 
ground). Furthermore, the static spacetime can be easily seen to be conformally flat and therefore it 
must necessarily be de Sitter (anti-de Sitter) |f[q| . Thus, the following result holds 

Proposition 1 Any static and axially symmetric spatially compact region in a de Sitter background 
must have the same geometry as the exterior manifold. 

In other words, de Sitter (anti-de Sitter) spacetimes do not admit non-trivial static and axially sym- 
metric spatially compact regions. Thus, we can assume from now on that that the FL spacetime is 
not de Sitter or anti-de Sitter, i.e. (p + p) FL ^ 0. 

If the static region is imposed to have an energy-momentum tensor with p + p r = (this includes 
the important subcases of vacuum, A-term and non-null Einstein-Maxwell solutions, see below), then 
( [42| ) immediately implies f — and r t = (<=> r(t) = r ). Thus, the matching hypersurface 
must be a two-sphere comoving with the cosmological flow. Conversely, if the matching two-sphere is 
comoving with the cosmological flow (/ = , r(t) — ro), then ( ff2| ) implies p + p r = 0, which shows 
the equivalence of the two conditions. In this situation, the differential equation for T tt (|3^) can be 
integrated to give 

T t = 



E 2 r - a 2 E 2 

(the integration constants are irrelevant, as discussed above) so that the line-element in the static 
region becomes 



^ 2 £ 2 ) dT 2 + °'* S 2 dt 2 + a 2 Y? (de 2 + sin 2 9d$ 2 ) 



S 2 r - a 2 E 2 

which is obviously spherically symmetric. Performing the coordinate change 

a(i)£(r ) = f 
and defining the standard mass function m(r) Jl9| as 

2 2m (a£) 



(45) 



jE 3 



the line element (W5h transforms into 



ds 2 =-(l- ^\ dT 2 + — ^_ + r 2 (dO 2 + sin 2 0# 2 ) . (46) 

f 

It is interesting to notice that this line-element is the most general static and spherically symmetric 
metric with p + p r = 0. Thus, all these metrics can be matched with an appropriate exterior FL 
spacetime through a spherically symmetric matching hypersurface. The energy-momentum contents 
and the only non-zero Weyl spinor coefficient, ^2, of this spacetime read 

2m f m ff m f f 2m f m 

p = -p r = PT = 2_ 5 ^ 2 = — ^_ + _^__ (47) 

ij*4d iy» xyy oT* 



while the energy-density and pressure in the FL cosmology are 

FL _ 6m FL _ 2m,r 

P ~~ ~^T' P 

Thus, we can conclude 
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Proposition 2 The matching between a static and axially symmetric spacetime and a Friedman- 
Lemaitre cosmology across a hypersurface described by a two-sphere comoving with the cosmological 
flow must have a spherically symmetric interior with metric given by where 771(f) is an arbitrary 
junction describing the total mass contained in the static region. 

Since vacuum and A-term spacetimes are particular cases of p + p r = 0, the two following results 
hold 

Corollary 1 The only A-term static and axially symmetric spatially compact region within a FL 
spacetime must be a two-sphere comoving with the cosmological flow and interior geometry given by 
with m — tiiq + (tuq and po constants). 



Corollary 2 The only static and axially symmetric, spatially compact, vacuum region in a Friedman- 
Lemaitre spacetime must be a two-sphere comoving with the cosmological flow and with Schwarzschild 
interior geometry (the Einstein- Straus model) 

Let us next discuss Einstein-Maxwell spacetimes, i.e. spacetimes with an energy-momentum tensor 
corresponding to an electromagnetic field outside any mass and charge sources. Static spacetimes with 
a pure null electromagnetic source are impossible |20|j . so we are left with the non-null case. In this 
case (and due to the diagonal form of the energy-momentum tensor) the so-called Rainich conditions 
(which are necessary and sufficient conditions for non-null electrovacuum spacetimes) fl8f , imply that 
(T a p) must take the form 

T™ = diag (77,-77,77, 77). 

where r\ is some positive quantity. Hence, p + p r = follows and the metric must take the spherically 

symmetric form (Efy). The general Einstein-Maxwell solution for such a metric is the well-known 

— 2 
Reissner-Nordstrom solution, with m(f) = M — (where M is the mass and e is electric charge of 

the point source). The scale factor and the energy contents in the FL spacetime are then given by 



2M e 1 



FL _ 6M 3e 2 



P 



a 3 £ 3 a 4 £ 4 



FL 



P 



a 4 £ 4 



Summarising, we have 



Proposition 3 The only static and axially symmetric electrovacuum region in a Friedman- Lemaitre 
spacetime must be a two-sphere comoving with the cosmological flow and with Reissner-Nordstrom 
interior geometry. 

Having discussed vacuum, A-term and electromagnetic fields, let us next analyse the perfect-fluid 
case. Imposing p r — px — and using the expressions above, we obtain a linear relation in cos 9. The 
vanishing of the coefficient in cos 9 implies necessarily / = while the other coefficient provides a 
rather long equation relating a(t) and r(t). Thus the following result also holds 

Proposition 4 Any static and axially symmetric perfect-fluid region immersed in a Friedman- 
Lemaitre cosmology must be spherically symmetric and its boundary must be an expanding (or con- 
tracting) two-sphere with its centre being at rest with respect to the cosmological flow. 

Here the terms expanding and contracting are defined with respect to the cosmological fluid flow. 

All these results above show that the most commonly used energy-momentum tensors imply nec- 
essarily that the matching hypersurface must be a two-sphere with its center being comoving with 
the cosmological fluid (and hence the matching hypersurface is spherically symmetric). Furthermore, 

1 o 

the interior geometry must also be exactly spherically symmetric in the relevant open set n (>V) 
(this set may be empty in some marginal cases - see the previous section- but even in that particular 
situation the conclusions on the form of the matching hypersurface hold). Consequently, spherical 
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symmetry is indeed an essential ingredient for the existence of reasonable static cavities in a FL ex- 
panding background. This implies that the Einstein-Straus model represents a very isolated situation, 
which poses serious doubts on the suitability of this model for providing a definitive answer to the 
problem of the influence of the cosmological expansion on the local physics. Even though this effect 
is probably very small at solar system scales (there is no experimental evidence for such an effect), 
this is not so clear at larger scales (ranging from galaxies to clusters of galaxies or even larger). In 
any case, more adequate models describing the problem should be analyzed before reaching definitive 
conclusions on this subject. 

Before finishing this paper, a final comment may be adequate. Solving any matching problem 
has always a dual interpretation. Due to the symmetry (+) «-» (— ) of the matching conditions, 
solving a matching between two spacetimes produces always as a side effect the matching of the 
two complementary regions discarded in the original matching (see |2l| l for a clear discussion of this 
fact). In our case, the complementary situation corresponds to an interior spatially homogeneous and 
isotropic perfect fluid (with FL metric) and an static, axially symmetric exterior. Thus, irrespectively 
of any conditions on the static exterior (for instance, no conditions on asymptotic flatness or similar), it 
turns out that the shape of the interior fluid must be a two-sphere at each instant of time. Furthermore, 
for the most relevant exterior metrics (in particular vacuum and electrovacuum) , the shape of the 
autogravitating fluid must be exactly spherically symmetric and the exterior geometry must be also 
exactly spherically symmetric. This side result is interesting in its own and can be phrased as follows. 

A Friedman-Lemaitre spacetime acting as a source for static and axially symmetric spacetimes with 
physically reasonable energy-momentum tensors (when the two canonical times match properly) re- 
quires that both the boundary of the fluid and the exterior geometry are exactly spherically symmetric, 
and hence the global spacetime must also be spherically symmetric. 
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